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Screening of Hydrodynamic Interactions in 
Semidilute Polymer Solutions: A Computer Simulation Study 

Patrick Ahlrichs, Ralf Everaers and Burkhard Diinweg 

Max— Planck— Institut fur Polymerforschung, Postfach 3148, D— 55021 Mainz, Germany 

We study single-chain motion in semidilute solutions of polymers of length A^ = 1000 with excluded- 
volume and hydrodynamic interactions by a novel algorithm. The crossover length of the transition from Zimm 
(short lengths and times) to Rouse dynamics (larger scales) is proportional to the static screening length. The 
crossover time is the corresponding Zimm time. Our data indicate Zimm behavior at large lengths but short 
times. There is no hydrodynamic screening until the chains feel constraints, after which they resist the flow: 
"Incomplete screening" occurs in the time domain. 

PACS Numbers: 83.10.Nn, 05.40.Jc, 47.11.+J, 61.25.Hq, 36.20.Ey 
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The dynamics of polymer chains in solution []1|,^ has been 
the subject of long-standing theoretical investigations, even 
for the simple case of flexible uncharged chains in good sol- 
vent. While in the dilute limit the validity of Zimm scaling 
predictions ^^ is generally accepted, as confirmed by ex- 
periments |y,^|^] and computer simulation [0-Q|, the theoret- 
ical understanding becomes much more involved as soon as 
one considers finite concentrations [p^|-p^. This is so due 
to the complicated interplay between excluded-volume inter- 
actions, hydrodynamic interactions, and entanglements. Only 
for the opposite limit of dense melts, where the first two inter- 
actions are fully screened, there exists a well-controlled de- 
scription in terms of the Rouse or reptation model |jl^,n8ltt9p . 
However, the details of the crossover, the underlying mecha- 
nism of the screening of hydrodynamic interactions, and the 
concentration dependence of the screening length have been a 
subject of considerable debate. 

In this Letter, we present the first computer simulation study 
which is able to contribute to the resolution of these ques- 
tions. Experiments, such as light scattering 120,2l|] or non- 
equilibrium methods [O] usually focus on collective concen- 
tration fluctuations, while single-chain motions are only ac- 
cessible by labeling techniques (neutron | jl4| ] or light [q]| scat- 
tering). Computer simulations can in principle analyze both 
types of motion; however, for reasons of statistical accuracy 
we had to confine ourselves to single-chain motion. 

We study the equilibrium fluctuations of a three- 
dimensional semidilute system of flexible bead-spring poly- 
mer chains with full excluded volume interactions, coupled 
to a hydrodynamic background to fully take into account hy- 
drodynamic interactions, using an efficient method which we 
have recently developed and tested [g]. The polymer system 
is simulated by Langevin stochastic dynamics, the solvent by 
a stochastic D3Q18 lattice Boltzmann model [E3pJ], and a 
point-particle coupling is introduced via a monomeric fric- 
tion coefficient. The present work uses the same model with 
the same parameters but in the semidilute regime. One par- 
ticular advantage, without which the study would have been 
unfeasible, is the fact that the lattice Boltzmann solvent does 
not alter the good solvent statistics of the chain conformations. 



We therefore first equilibrated the multi-chain system without 
the computationally expensive solvent, using a combination of 
stochastic dynamics and slithering-snake Monte Carlo moves 
comprising several momomers (roughly one blob, see below). 
This run produced a set of configurations, which were after- 
wards coupled to the solvent. For further details of the model, 
we refer the reader to Ref. |Q]. 

Semidilute systems are characterized by a very low 
monomer concentration c, which is nevertheless large enough 
to induce strong overlap of the coils. The static conforma- 
tions are well understood [Q] in terms of the "blob size" <^5, 
i. e. the typical mesh size of the temporary network. On 
scales below ^5, the chains are self-avoiding walks (SAWs) 
characterized by the scaling law R ^ aN^ , where a is the 
monomer size, R the chain extension, N the degree of poly- 
merization and V « 0.59. The concentration dependence of 
^5 results from c '^ S,s'^{is/o)^^" ■ On scales above 1^5 the 
excluded-volume interaction is screened, and the chains are 
random walks (RWs, R ^ aN^^"^). The overall chain is 
thus a RW of blobs with R^ - (s^/i^s/aY'" ■ This pic- 
ture implies that rather long chains are necessary in order to 
clearly observe both regimes; guided by the idea of having 
roughly 30 blobs of 30 monomers each available, we chose 
N = 1000, and varied ^s by studying the concentration val- 
ues c == 0.00837, 0.0397, 0.0734, 0.134, 0.231 for the statics, 
and the latter three values for the dynamics (data are always 
given in the Lennard-Jones unit system of Ref. [^). The num- 
ber of chains M = 50 was kept fixed; this is large enough to 
ensure that even the most concentrated system does not exhibit 
self-overlap due to the periodic boundary conditions. Our 
data for the static chain conformations are in perfect agree- 
ment with the blob scaling picture p5[|, as have been those of 
previous extensive Monte Carlo simulations pdETn. 

Dynamic scaling for a single chain which exhibits no spe- 
cial length scale except a and R implies t (x R^, where r is 
the conformational relaxation time, and z — 3 for the Zimm 
model (applicable to dilute solutions without chain overlap 
where hydrodynamic interactions are fully developed), while 
z ~ 4 for the RW Rouse model (applicable to dense melts 
where hydrodynamic interactions are fully screened; we do 
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FIG. 1. Scaled /Ar^\ for different concentrations. 



not consider reptation-like slowing down, which occurs only 
for longer chains, and does not play any role for our simu- 
lation data). In terms of the chain diffusion constant D this 
corresponds, via Dt ^ R^,to D (x R^^ for Zimm (chain as 
a Stokes sphere) and D ex N^^ for Rouse (monomers as inde- 
pendent friction centers). Furthermore, the scaling of lengths 
with the corresponding times implies a subdiffusive behavior 
for the single-monomer mean square displacement (A?-'^) ex 
i^/^, and a k'^t behavior for the single-chain dynamic struc- 
ture factor 5(fc, t) = A^-^ ^^^. /exp (ik ■ {fi{t) - fj{0))] \ 

in the scaling regime of intermediate length scales (between a 
and R) and time scales (between tq, the microscopic time for 
monomer relaxation, and r). 

For a semidilute system, one expects a crossover between 
these cases. Indeed, our data for ( Ar^) do exhibit a crossover 
from a Zimm-like t^^^ behavior at short times to t^^^ at 
longer times. The behavior at short length and time scales 
is thus Zimm-like. The pure Zimm model [|l|] predicts that 
the decay rate, i. e., in the given context, the prefactor A 
of the law (Ar^) ~ At^/^ should only depend on solvent 
viscosity and not on concentration. We nevertheless found a 
weak concentration dependence of A (roughly 20% within the 
given concentration range, see below). Figure |]thus studies 
(Ar^) /{At^^^) = f{t/tc), where tc is the concentration- 
dependent crossover time, which again is the Zimm relax- 
ation time of a dynamic crossover length (or hydrodynamic 
screening length), tc oc ^|^, and f{x) oc x"^/^ for large x. 
We find a very good data collapse assuming that £^h oc ^s 
or tc oc c^^''/'^"^^' = c~^-^, as done in Fig, p]. The as- 
sumptions ^// oc c~^ [ pXipJi l and (^h oc c~^/^ [|16[| produced 
significantly poorer collapses. The simulation thus confirms 
de Gennes' prediction ^h oc S,s [0^ as also observed in ex- 
periments (except for weak corrections pO| , p2| |). 

S{k, t) was studied for the most dilute system c = 0.0734. 
We prefer scaling plots of the raw data using asymptotic expo- 
nents over fits to functional forms derived from approximate 
theories. From S{k, 0) we estimated the crossover wavenum- 
ber as kc w 0.45 and the scaling regime as 0.2 < k < 1.5. 
From (Ar^) we read off tc ~ lO'^; the nonuniversal regime 
t < 30 was discarded. Figures and pi show that for short 
times t <^tc the decay can be described quite well by Zimm 
scaling, regardless of wavenumber, while for t ^ tc there 
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FIG. 2. Decay of S{k,t)/S{k,t) on length and time scales as 
indicated, suggesting Zimm scaling for short times (main figure) and 
Rouse scaling for long wavelengths (inset). 



is a simultaneous smooth crossover to Rouse dynamics for 
those wavenumbers which have not yet fully decayed, i. e. for 
k < kc- Note that the initial Zimm regime of these wavenum- 
bers can be easily overlooked in the representation of the inset 
of Fig. I 

The physical picture which results from this observation is 
thus free Zimm motion up to the crossover time (on all length 
scales), after which screening sets in, leading to Rouse-like 
motion. Hence the most important finding of our simulation is 
that hydrodynamic screening must necessarily be viewed as a 
dynamic time-dependent phenomenon. We consider this to be 
the logical consequence of the (correct) original treatment by 
de Gennes [ ]ll| ] (see below). Nevertheless, this has so far been 
overlooked in the literature, the main reason being that single- 
chain motion on length scales beyond ^ is not accessible to 
standard scattering experiments [ |2C| , pi| ] which are sensitive to 
collective concentration fluctuations: On scales fc<^ <C 1 the 
overall solution is homogeneous, and one observes a simple 
diffusive decay cxp{—Dcoopk^t) with Dcoop oc ^~^. Accord- 
ingly, single-chain motion on scales beyond £, was not treated 
explicitly in Ref . 1 1 1 1 . The experiments on labeled chains 
p,n4| produced data which are fully consistent with our view, 
but were interpreted incorrectly (see below). The rest of the 
paper will be devoted to further discussion, and elucidation of 
the underlying mechanism. 

Hydrodynamic interaction is the presence of long-ranged 
correlations in the stochastic displacements of a system of 
Brownian particles, caused by fast diffusive momentum trans- 
port through the surrounding fluid. They can be calculated 
by solving the stationary Stokes equation around a system of 
spheres p^], resulting in a complicated multipole expansion 
which contains many-body terms representing the multiple 
scattering of the flow. In the dilute limit, however, it is suffi- 
cient to just consider the leading-order pair interaction, which 
decays like r^^, where r is the interparticle distance (Oseen 
tensor). Conversely, screened hydrodynamic interactions are 
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FIG. 3. S{k, t)/S{k, 0) in the RW regime 0.2 < fc < 0.45 witli 
both Zimm and Rouse scahng, using a representation which empha- 
sizes the short-time behavior. 



described by a Yukawa-like decay r~^ exp(— r/^/f ) defining 
the hydrodynamic screening length S^h- Such an interaction 
occurs for Darcy flow through a porous medium, where fixed 
frictional obstacles with friction constant ( exert a force —(u 
on the flow with velocity u. Denoting the obstacle concen- 
tration with c, the flow is described, on scales beyond the 
typical interparticle distance, by a modified Stokes equation 
pdu/dt — rjW^u — C,cu, where p is the fluid density and 77 the 
viscosity. This implies rj^^" = (c. 

The simplest approach to hydrodynamic screening in poly- 
mer solutions just replaces the r^^ Oseen interaction by a 
screened Yukawa-like interaction, leading to uncorrected dis- 
placements of monomers whose distance exceeds £^h- The 
resulting motion of the chain is Zimm-like on short length 
and time scales and Rouse-like on length scales beyond ^h 
for all times [|l4|]. The Darcy flow thus produces the desired 
crossover. 

Unfortunately, this picture generates as many questions as 
it answers. In particular, the obstacles must be the mobile 
polymer chains themselves, whereas strict Darcy flow requires 
fixed obstacles. Moreover, momentum is present and being 
transported infinitely far in polymeric as well as in simple flu- 
ids. This fundamental conceptual difficulty was recognized 
by Richter et al. [[l4|]. In their "incomplete screening" model 
they proposed that the hydrodynamic interaction should cross 
over to a second r~^ regime on very large scales, but with 
the overall viscosity as a prefactor For the single-chain 
short-time behavior this model also predicts Rouse-like mo- 
tion. However, this regime is now restricted to length scales 
^H < k^^ < S,HVsoiution/risoivent- On larger scalcs there is 
an additional Zimm regime. Richter et al. [ |l4| ] used this to in- 
terpret the mixture of Rouse- and Zimm-like signals in their 
scattering data. Similar arguments were used by Martin 0], 
who observed Zimm scaling on all length scales in the initial 
decay rate of dynamic light scattering of labeled chains. 

The simple model and the more refined version by Richter 
et al. [pp are at variance with both our data and our theoretical 



arguments for the short-time behavior, see below. It should be 
noted that the "incomplete screening" model must have fun- 
damental conceptual difficulties, since rjsoiuUon appears in the 
short-time dynamics, although rjsoiuUon is established only 
on time scales beyond the overall chain relaxation. 

It is therefore clear that a consistent theoretical description 
has to study the dynamics of the coupled polymer-solvent sys- 
tem. The first attempt by Freed and Edwards [HQ] considered 
a multiple scattering series of the flow around the monomers, 
which is in spirit somewhat similar to the multipole expan- 
sion of Ref. [|2^]. After some approximations an effective 
Darcy equation results, with ^h oc c^^. Later this scheme 
was shown to be inadequate; in Ref. [|l2|| it was proven that 
a system of phantom chains (which do not interact with each 
other, but to which the original approach [ |l(]| ] should apply 
as well) does not exhibit any hydrodynamic screening what- 
soever. This is consistent with new results for colloidal sus- 
pensions, where possible hydrodynamic screening has been 
discussed recently; the result of large-scale computer simula- 
tions [ p9| ] was the absence of screening. 

With respect to hydrodynamic screening in polymer solu- 
tions we can thus draw the following conclusions: (i) nei- 
ther does the presence of higher-order terms of the multipole 
expansion [|8| at finite concentrations of scattering centers 
lead to screening; (ii) nor can such terms be of any impor- 
tance in the semidilute limit, where one can reach arbitrar- 
ily small monomer concentrations c, while still keeping the 
chains at strong overlap; (iii) as this "colloidal" mechanism of 
screening does not apply, the underlying physics must rather 
be polymer-specific; and (iv) the mechanism must lead to a 
time-delayed screening. 

Concerning the short-time behavior, we note that the 
semidilute system is governed by a Kirkwood diffusion equa- 
tion [11]], with Si pure Oseen-type r^^ diffusion tensor, which 
describes the short-time diffusive behavior, and a force term 
due to connectivity, excluded volume, and entanglements. 
Within this formalism, it can be shown rigorously [|l|] that 
the initial decay rates of correlation functions are only gov- 
erned by the diffusion tensor and the statistics of the chain 
conformations. In particular, considering the initial decay rate 
of the single-chain structure factor [|[], one obtains the same 
formula as for an isolated chain in solvent — the effect of 
the other chains is merely the modification of the conforma- 
tions. Zimm chains, however, have always z ~ 2> indepen- 
dently of chain statistics; the initial decay rate is given by 
V{k) ^ {kBT/ri)k^, while the fractal dimension only en- 
ters the prefactor [nl] . For systems in the semidilute limit we 
thus conclude, in accordance with our simulation results and 
the experimental data by Martin |a], that for short times the 
single-chain dynamics is Zimm-like, independently of length 
scales. 

In his pioneering 1976 paper [ pT[ | de Gennes noticed that 
the decisive mechanism for screening is the connectivity and 
the strong coupling to the temporary matrix (expressed in 
terms of an elastic gel, which is physically more appropri- 
ate than a rigid porous medium): After the time needed for 



a blob to move its own size, which is the blob's Zimm time 
tc, the chain will, on average, feel the constraints by the tem- 
porary matrix. From then on it is unable to follow the flow, 
but rather lags behind, and starts to produce Darcy-type fric- 
tional resistance. As the blob can be envisioned as a Stokes 
sphere with radius ^ and friction coefficient ^ -q^, the obsta- 
cles which produce the Darcy flow are not the monomers but 
rather the blobs. Hence, the hydrodynamic screening length 



is given by 77^^ 
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77^5^0 , i. e. the hydro- 



dynamic screening length is, apart from prefactors, identical 
to the static screening length, S^h oc Cs- This argument [ pl| 
makes the picture fully self-consistent. On length scales be- 
yond ^, and time scales beyond tc, the semidilute solution is 
just a Rouse melt of blobs, while the conformations within 
the blobs are already fully relaxed. In the Rouse regime, mo- 
mentum transport is no longer described by a simple Navier- 
Stokes equation. It rather occurs mainly along the chain back- 
bones, due to the connectivity forces. This results in a very 
efficient randomization of a locally applied "kick". 

It should be noted that our simulated system deviates some- 
what from that ideal scenario. The most dilute system has a 
density of 9% of a typical dense melt, and thus one must ex- 
pect that higher-order terms in the multipole expansion [ p8[ | 
do play a role. We believe that these are the main source of the 
c dependence of the prefactor of the initial t^/^ law in (Ar^). 
We expect finite size effects [^g] with respect to the linear box 
size L to be rather small, since for our data (kL)^^ < 0.06. 

To summarize, we have presented the first computer sim- 
ulation study which was able to study the dynamic crossover 
from Zimm to Rouse behavior in semidilute polymer solu- 
tions. This was made possible by a novel algorithm, whose 
essential feature is the replacement of the solvent by a Navier- 
Stokes background which is coupled dissipatively to the 
monomers. Our results are fully consistent with de Gennes' 
scaling picture [|ll|], and emphasize the fact that hydrody- 
namic screening is a dynamic effect which becomes relevant 
only after the crossover time. Incomplete screening thus in- 
deed occurs, however not on large length scales [Q, but on 
short time scales. Any theoretical description which builds 
upon a screened hydrodynamic interaction depending only on 
distance, but disregards the time dependence, cannot describe 
the phenomena correctly. 

We gratefully acknowledge fruitful discussions with K. 
Kremer, and thank the Max Planck Society for generous al- 
location of Cray T3E computer time at RZG. 
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